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ABSTRACT 

We compute the X-ray spectra produced by non-static coronae atop accretion discs around 
black holes and neutron stars. The hot corona is radiatively coupled to the underlying disc 
(the reflector) and generates an X-ray spectrum which is sensitive to the bulk velocity of the 
coronal plasma, (3 — v/c. We show that an outflowing corona reproduces the hard-state spec- 
trum of Cyg X-1 and similar objects. The dynamic model predicts a correlation between the 
observed amplitude of reflection R and the X-ray spectrum slope F since both strongly depend 
on (3. A similar correlation was observed and its shape is well fitted by the dynamic model. 
The scattering of soft radiation in an outflowing corona can also account for the observed 
optical-UV polarization pattern in active galactic nuclei. 

Key words: accretion, accretion discs - radiative transfer - gamma-rays: theory - galaxies: 
Seyfert - X-rays: general - stars: individual: Cyg X-1 



1 INTRODUCTION 

The hard X-ray spectra of galactic black holes (GBHs) and active 
galactic nuclei (AGN) indicate the presence of hot plasmas with 
temperature kTe ~ 100 keV and scattering optical depth tt ~ 1 
in the vicinity of accreting black holes (see reviews by Zdziarski 
et al. 1997; Poutanen 1998). The plasma can be identified with 
a corona of a black hole accretion disc (e.g. Bisnovatyi-Kogan 
& Blinnikov 1977; Liang 1979; Galeev, Rosner & Vaiana 1979; 
see Beloborodov 1999a for a recent review, hereafter B99a). The 
corona is likely to form as a result of magnetorotational instabili- 
ties in the disc and the buoyancy of the generated magnetic field 
(Tout & Pringle 1992; Miller & Stone 2000). The corona is prob- 
ably heated in flare-like events of magnetic dissipation producing 
the variable X-ray emission. 

The dominant cooling mechanism of the flaring plasma is 
Compton cooling. The observed hard X-rays are generated by the 
Comptonization process, i.e. by multiple upscattering of seed soft 
photons by the hot electrons in the corona. The Comptonization 
generally produces power-law X-ray spectra. In addition to the di- 
rect power-law radiation from the corona, one also observes the 
X-rays reflected by the underlying (relatively cold) accretion disc 
(e.g. George & Fabian 1991). The arising reflection features in the 
spectrum, in particular, the Compton bump and the fluorescent iron 
line provide diagnostics for accretion models. 

The cold accretion disc partly reemits the incident X-rays in 
the form of soft thermal radiation. This radiation cools the corona. 
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providing the feedback loop which regulates the temperature of the 
corona (Haardt & Maraschi 1993, hereafter HM93). The geometry 
of the corona can hardly be derived from first principles. It might 
be a large cloud covering the whole inner region of the disc. It may 
also be a number of small-scale blobs with short life-times, with the 
energy release concentrated in space and time. The resulting X-ray 
spectrum is however not sensitive to the exact shape of the cloud, 
its density distribution, and other details. The only important pa- 
rameter is the effective feedback factor (HM93; Stern et al. 1995b; 
PS96; Svensson 1996) that is the fraction of the X-ray luminosity 
which reenters the source after reprocessing. 

Previous computations of the disc-corona models all assumed 
that the corona is static (e.g. HM93; Stem et al. 1995b; Poutanen 
& Svensson 1996, hereafter PS96; see Svensson 1996, Poutanen 
1998 for a review). The model was successfully applied to Seyfert 1 
AGN, however, it was found to disagree with observations of some 
black-hole sources in the hard state, for instance, Cyg X-1 (see e.g. 
Gierlinski et al. 1997 and Section^). Three alternatives have been 
suggested: (i) the cold disc is disrupted in the inner region (e.g. 
Poutanen, Krolik & Ryde 1997; Esin et al. 1998), (ii) the disc is 
highly ionised (Ross, Fabian & Young 1999; Nayakshin, Kazanas 
& Kallman 2000), and (iii) the coronal plasma is moving away from 
the disc and emits beamed X-rays (Beloborodov 1999b, hereafter 
B99b). We here study the third scenario. 

B99b showed that the static model is not self-consistent for 
-dominated flares and argued that the hot plasma should be 
ejected with a mildly relativistic bulk velocity l3 = v/c ^ 0.5. 
The bulk acceleration is also efficient in flares dominated by nor- 
mal e — p plasma if the compactness parameter of the flare exceeds 
^ 100 (B99a). The likely pattern of ejection resembles that of solar 
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flares (scaled to higher velocities because of higher compactness). 
The ejection velocity may be directed away or towards the disc (and 
it may change) with the preferential direction away from the disc. 

Mildly relativistic bulk motion causes aberration of the X-ray 
emission and strongly affects both the amplitude of reflection R 
and the spectrum slope V. Using a simple analytical model, B99b 
estimated the dependence of R and F on /3 and found that the 
low amplitude of reflection _R ^ 0.3 and the spectrum hardness 
r ~ 1.6 observed in the hard state of Cyg X-1 can be explained 
with 13 ~ 0.3. The model also explains the observed correlation 
between R and V (B99a; Zdziarski, Lubihski & Smith 1999; Gil- 
fanov, Churazov & Revnivtsev 2000). 

In the present paper, we perform exact computations of the 
X-ray spectra produced by dynamic coronae. We use a non-linear 
Monte-Carlo code (Malzac & Jourdain 2000) which is based on the 
large-particle method (Stern 1985; Stem et al. 1995a). The formu- 
lation of the problem and details of simulations are described in 
Section ^. In Section ^ we compute the spectra in the static case 
which was well studied previously and illuminate the problem of 
static models. In Section ^ we present spectra from dynamic coro- 
nae. The results are compared with observations in Section H. 




Figure 1. A heated cylinder atop the accretion disc. The cylinder is divided 
into nine computational cells of equal volumes. The hot plasma moves 
through the cylinder with a velocity (3. 



2 SETUP 

2.1 A heated cylinder 

Consider a cylinder of radius r and height h located atop the accre- 
tion disc (Fig. 1^). The cylinder may be associated with a hot outflow 
covering the disc or a heated magnetic tube in a compact flare. In 
the limit r S> /i we get a slab geometry of the corona, which has 
the largest feedback. In the opposite case, r <gi h, the feedback 
tends to zero (most of the reprocessed radiation goes away from 
the cylinder). We study the sequence of geometries parametrized 
with the ratio h/r. The plasma in the cylinder is assumed to have 
a constant density and it is heated homogeneously with a constant 
rate. 

The plasma moves through the cylinder with a velocity (5 di- 
rected normally to the disc. We assume that the heated cylindri- 
cal volume is static with respect to the disc. The plasma motion 
then implies a flux of particles through the cylinder. If /3 < then 
the particle flux is absorbed by the disc at the bottom. If /3 > 
then the particles escape through the top of the cylinder. Outside 
the heated region, the particles are immediately (on a time-scale 
<C /i/c) cooled down to the Compton temperature of the radiation 
field, kTc ~ 1 — 10 keV. In our simulations, we do not take into 
account the scattering on the escaped cold particles. Their density 
is likely to be reduced if the plasma flows out along diverging mag- 
netic lines. 

In this paper, we restrict our consideration to thermal coronae 
and assume that the heated electrons have a Maxwellian distribu- 
tion with a temperature Te in the plasma rest frame. The cylinder 
is divided into nine cells of equal volumes (see Fig. |l]) and the tem- 
perature Te is calculated from the local heating=cooling balance 
in each of the cells separately. Note that the resulting equilibrium 
temperature is not homogeneous. 



of the coronal plasma and neglect soft radiation generated vis- 
cously inside the accretion disc. This should be a good approxi- 
mation if the corona releases energy in the form of strong concen- 
trated flares; then, locally, the soft flux from inside the disc is much 
smaller than the reprocessed fluxQ We also neglect additional pho- 
ton sources such as bremsstrahlung and cyclo-synchrotron emis- 
sion in the corona (see e.g. Wardzihski & Zdziarski 2000). An in- 
dependent argument in favour of the reprocessed radiation as the 
main cooler comes from the observed R — T correlation (Zdziarski 
et al. 1999; Gilfanov, Churazov & Revnivtsev 1999; see Section|5|). 

We assume that the reflecting material of the disc is suffi- 
ciently dense so that the ionisation parameter ^ 10^ and the 
ionisation effects are weak (e.g. Zycki et al. 1994). Then the albedo 
is small, a ~ 0.2, and most of the X-rays impinging the disc 
are reprocessed. We assume that the reprocessed flux has a quasi- 
blackbody spectrum (possibly diluted) with a constant temperature 
Tbb. In the simulations, we consider two cases: fcTbb = 150 eV 
and fcTbb = 5 eV, representing the typical blackbody temperatures 
in GBHs and AGN, respectively. 

The feedback factor D is defined as the ratio of the soft black- 
body luminosity entering the cylinder, Ls, to the total luminosity 
outgoing from the cylinder, L. The Comptonization time-scale is 
quite short, ~ {h/c) log(meC^/3fcTbb), and we assume the 
heating rate to be steady on time-scales t ^ t,. Then the X-rays 
are generated in a quasi-stationary regime and the feedback fac- 
tor determines the Compton amplification factor of the hot plasma 
A — L/Ls through the relation 



DA= 1. 



(1) 



The amplification factor in its turn is the main magnitude control- 
ling the spectral slope F of the Comptonized radiation (see B99a 
and Section 4.3). 



2.2 Feedback 

The X-rays from the corona strike the underlying disc and get re- 
processed. We assume that reprocessed radiation is the main cooler 



* Yet, the total surface-integrated intrinsic luminosity of the cold disc does 
not need to be much smaller than the coronal luminosity since the covering 
factor of the corona may be small (Haardt, Maraschi & Ghisellini 1994). 
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2.3 Radiative transfer and energy balance 

We compute the plasma temperature simultaneously with the ra- 
diative transfer. Our code is based on the non-linear Large Particle 
Monte-Carlo method (Stern 1985; Stern et al. 1995a). In contrast 
to the standard Monte-Carlo technique, this method allows one to 
follow the path and successive interactions of photons and particles 
in parallel. The code is described and tested in Malzac & Jourdain 
(2000). 

We start a simulation from an initial (non-equilibrium) state 
and follow the evolution of the system plasma plus radiation un- 
til a steady state is achieved. The time step is about 10^'^^/c. We 
compute the evolution of the temperature distribution in the cylin- 
der from the local energy balance. The plasma gains energy owing 
to constant heating and loses energy via Compton cooling. The dif- 
ference between heating and cooling determines the change of the 
temperature during each time step. 

The radiative transfer (photon tracking) is dealt in the lab 
frame in the same way as for static coronae (Malzac &. Jourdain 
2000). The only difference is that the velocity of the scattering 
electron is generated from an isotropic Maxwellian distribution in 
the plasma rest frame and then it is Lorentz transformed to the lab 
frame. 

In all simulations in this paper, we assume a constant density 
throughout the cylinder and parametrize our models by Thomson 
optical depth tt = n^GTh rather than a dissipation rate Ldiss 
or compactness /diss = {Ldiss / h){(JT /me(?). Strictly speaking, 
the model is self-consistent for low temperatures or low compact- 
nesses, when pair production is not important. In pair dominated 
blobs, the density cannot be assumed homogeneous: the density 
distribution is then determined by the local pair balance. Yet, mod- 
els with pair production give similar results if they have same tt. 

The model has four parameters: (i) Thomson optical depth tt 
(defined along the height of the cylinder), (ii) height to radius ra- 
tio of the cylinder h/r, (iii) bulk velocity and (iv) blackbody 
temperature Tbb. Given these parameters the code computes the 
emitted spectrum as a function of the inclination angle i (the angle 
between the disc normal and the line of sight). We consider 10 an- 
gular bins of equal width Acqs i = 0.1. The spectra are averaged 
within each bin. 



2.4 The intrinsic and reflected components 

For each inclination we compute two components of the observed 
emission: 

(i) The 'intrinsic' Comptonised X-rays, i.e. the photons coming 
to the observer directly after the upscattering in the hot blob. The 
slope of the intrinsic X-ray spectrum (the photon index F) is evalu- 
ated from a least- square fit in the 2-10 keV range. 

(ii) The reflected/reprocessed radiation from the underlying 
disc. We count only those reflected photons that are not scattered 
in the blob, thus accounting for the attenuation of reflection by the 
blob. 

We model the reflector as an infinite slab and compute the re- 
flection component assuming neutral reflecting material with stan- 
dard abundances (Anders & Ebihara 1982). The amplitude of re- 
flection is defined as the ratio of an observed reflected component to 
that expected from an isotropic point source illuminating the slab. 



The simulations yield the intrinsic Comptonized spectrum emit- 
ted at an inclination i and the reflected luminosity Li-ofl(*)- We 
then compute the reflected luminosity L|.^g(i) for an isotropic point 
source with the same intrinsic spectrum and find R{i) from equa- 
tion We fitted some of our simulated spectra with the PEXRAV 
model (Magdziarz & Zdziarski 1995) under XSPEC. The best fit 
values we found for R are very close to those derived using equa- 
tion Bt,. 



3 STATIC CORONAE 

We have computed a set of static models (/3 = 0) with differ- 
ent h/r, Tt, and Tbb, and determined R and F in the calculated 
spectra. The results are shown in Fig. ^. For simplicity, the pre- 
sented amplitude of reflection is averaged over inclinations, R = 

R{i)d{cosi). 



3.1 Feedbaclt, amplification, and the spectral index 

In the static case, the feedback factor -Dstatic is determined mainly 
by the geometry of the cloud, i.e. by h/r in our parametrization. 
We express the feedback factor as 

I^static = ^s',Lc = - - Ms), (3) 

where a is the energy-integrated albedo of the disc. This equation 
defines the effective geometrical parameter fi^ for a static cloud 
(B99a,b). The difference 1 — /^s is the fraction of reprocessed ra- 
diation that returns to the hot cloud. To the first approximation, 
can be associated with fio = /i/(4r^ + h'^Y''^, where 1 — /lo cor- 
responds to the angular size of the cylinder base measured from its 
centre. Numerically, we find fis ~ (4/5)/io (see Fig. |3|b). 

As one can see from Fig. |3| the effective /is is not fully deter- 
mined hy h/r but also depends on tt. The strongest dependence 
on Tt is observed for flat geometries (h/r ^ 1): fis and A sub- 
stantially decrease as tt increases. This effect can be understood 
analytically. Consider a slab corona and denote the radiation fluxes 
outgoing from the slab at the bottom and the top as F- and F+, 
respectively. The reprocessed flux entering the slab from below is 
Fs = {1 — a)F-. The amplification factor is, by definition. 



A = 



F-+F+ 1 + F+/F_ 



Fs 



(4) 



One can estimate the ratio F+ /F- by solving the radiative transfer 
equation in the Eddington approximation (see Rybicki & Lightman 
1979, ch. 1) and assuming isotropic scattering. Taking into account 
that the energy sources are homogeneously distributed in the slab, 
we get 



F- 



1 + 2/V3 
1 + VStt/2' 



Ms 



4(1 + VS) + Stt ' 



(5) 



R{i) = 



cfl(«) 



(2) 



These expressions are in good agreement with the results of simu- 
lations (see Fig. I h/r = 0). If Tt is large, the radiation is trapped 
at the bottom, near the reflector. In the limiting case tt ^ 1, one 
has F+ < F-, A ^ {1 - a)~^ and ^s -1. Then almost aU 
the X-rays are reprocessed because the source emits most of them 
towards the reflector. The bottom layers of the slab are cold and 
emit a soft spectrum for which the disc albedo a is close to zero 
(see Fig. ^) and A ^ 1. Such a small A results in the high F (see 
Fig. |c). 
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Figure 2. Characteristics of tlie static disc-corona model (/3 = 0). The active blob is a cylinder of height h and radius r. The model is determined by three 
parameters: h/r, the optical depth of the cylinder, tt (measured in the vertical direction), and the temperature of the soft radiation entering the blob, T^b- 
Filled and open symbols correspond to fcT^b = 150 eV (GBHs) and fcT^b = 5 eV (AGN), respectively. Squares and triangles show the cases tt = 0.5 and 
TT = 3. (a) The average electron temperature in the blob versus h/r. (b) The inclination-averaged amplitude of reflection R versus h/r. The dashed curves 
represent the analytical formula (^. (c) The inclination-averaged photon index of the intrinsic X-ray spectrum F (measured in the 2 — 10 keV band) versus 
h/r. (d) R versus F. 



At small Tt, the emitted spectrum hardens (Fig. gc). This is 
caused by the fact that the number of photons scattered in an opti- 
cally thin blob is small. When the Comptonized spectrum is made 
of fewer photons, the same luminosity can be emitted only if the av- 
erage photon energy is high. This leads to a hard spectrum (Fig. ^) 
with a break at high energies. The optically thin models have high 
temperatures (Fig. Ha): at fixed geometry and Tbb the product ttTe 
stays approximately constant (HM93; Stern et al. 1995b; Svensson 
1996; see also Fig. |l^ and hence an optically thin blob must be 
hotter. 



Blobs with high h/r ('elongated cylinders') have low feed- 
back and the resulting spectrum is hard (Fig. ^). This behaviour is 
known from previous simulations (see Svensson 1996 for a review; 
the previously computed spectra of detached spheres are similar 
to those of elongated cylinders). The product ttTe increases with 
h/r; for fixed tt it implies a growth of the plasma temperature 
with h/r (Fig.|[t). 

The emitted spectrum depends on Tbb- For this reason, we 
perform computations for GBHs (fcTbb = 150 eV) and AGN 
(fcTbb ~ 5 eV) separately. For the same h/r and tt, we get a 



© 0000 RAS, MNRAS 000, 000-000 



X-ray spectra of dynamic coronae 5 



10 



(a) 



■ ^ A 

■ □ A S 

■ □ . s 

■ □ . s 

■ □ ^ ^ 

■ ° , ^ 
. □ ^ 

■ ° A 



b5 











2 

h/r 



Figure 3. (a) Amplification factor, A, and (b) tlie effective geometrical factor (see eq. |3|) for static coronae. The symbols have same meaning as in Fig. 
The dashed curve shows (4/5)/io, the /iq is given in equation (n). 
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harder spectrum for GBHs than for AGN (despite the fact that the 
feedback factor and A are not sensitive to Tbb). The hardening is 
caused by two reasons: (i) The energy range in which the Comp- 
tonized luminosity is liberated is narrower in GBHs. Therefore, in 
order to radiate a given amount of luminosity, the spectrum has to 
harden and cut off at higher energy, (ii) Low-tt models are hot 
and have an anisotropy break in the spectrum (Stern et al. 1995b). 
In GBHs, this break is in the 2-10 keV band and affects the mea- 
sured r. 



3.2 The amplitude of reflection 

The reflection component produced beneath the blob is partly scat- 
tered and attenuated when passing through the blob. To the first 
approximation, the attenuation effect can be described analytically, 

7? = Mo + (l-Mo)e-^^'^-™\ Mo= , =■ (6) 

+ (V2)2 

Here it is assumed that the fraction /lo of the reflected luminosity 
goes directly to the observer. The remaining fraction, 1 — ^o, is 
reflected at the base of the cylinder and attenuated when transmitted 
through the hot plasma. The cylinder optical depth seen by soft 
photons emitted at the base depends on tt and h/r; it equals ~ tt 
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Figure 5. The effect of bulk motion on the emitted spectra. Here h/r = 2, tx = 3, and nearly face-on inclination is assumed, . 9 < cos i < 1 . The spectra 
are normalised to the blackbody peak. Left panel shows the case of AGN (fcTbb = 5 eV) and right panel corresponds to GBHs (fcTbb = 150 eV). 



for h/r = and vanishes at h/r ^ 1. We interpolate between the 
two extremes as rT(l — /io) and put this effective optical depth into 
the exponential in equation (^. 

The attenuation formula assumes that the scattered reflec- 
tion component is completely destroyed. Indeed, the 'Comptonised 
reflection' component has a power-law spectrum where the reflec- 
tion features (the Fe lines and the edge) disappear. 

Equation is in good agreement with the simulation results 
(see Fig. ^). The attenuation effect is strong if h/r is small (so 
that a large fraction of the reflected radiation has to pass through 
the cylinder) and if tt is large. In the optically thick regime, tt ^ 
3, almost all reflected photons reentering the cloud are Compton 
scattered and the attenuation effect saturates: further increase in tt 
does not reduce R. 

For elongated cylinders, the attenuation affects only a small 
fraction of the reflected luminosity and R approaches unity at high 
h/r. The reflection amplitude is higher in GBHs, especially in the 
high-temperature models (with small tt and large h/r) where R 
substantially exceeds unity (see Fig. ^). The explanation is as fol- 
lows. The radiation scattered once is strongly anisotropic and col- 
limated towards the reflector (HM93). In GBHs, because of rela- 
tively large temperature of seed soft photons, the first scattering 
order peaks above I keV and contributes significantly to the re- 
flected flux. As a result, R increases. By contrast, in AGN, the 
Comptonized emission above 1 keV is made of multiply scattered 
photons and the anisotropy effects are small. 



3.3 R — T anticorrelation 

Models with small h/r have soft spectra because of large feed- 
back. These models also have lower reflection amplitudes because 
of stronger attenuation (see Section 3.2). Therefore, variations in 
h/r produce an anticorrelation between R and F (Fig.|^). The an- 
ticorrelation disagrees with the data (see Section pb. A static model 
with a neutral reflector never predicts a small R simultaneously 
with a hard spectrum. In particular, the parameters R ~ 0.3 and 



r ~ 1.6 observed in Cyg X-I and similar hard-state objects cannot 
be reproduced. 



4 DYNAMIC CORONAE 
4.1 Example spectra 

In the dynamic case, /3 7^ 0, the corona spectrum is strongly af- 
fected by the bulk motion of the hot plasma. This is illustrated in 
Fig. ^ where we compare the cases /3 = 0.3 and f) = —0.2 with 
the static case. The spectra are computed for the same optical depth, 
Tt = 3, and geometry, h/r = 2. 

In the case of /3 = 0.3 (plasma moves away from the disc), 
the observed Comptonized luminosity is enhanced as a result of 
relativistic aberration. The X-rays are beamed away from the disc, 
and the reprocessed and reflected luminosities are reduced. The low 
feedback leads to a hard intrinsic spectrum. 

In the case of /9 = —0.2 (plasma moves towards the disc), the 
Comptonized luminosity is beamed towards the disc and the repro- 
cessed and reflected components are enhanced. The high feedback 
leads to a soft intrinsic spectrum. 

Since tt is fixed in Fig. ^, a high (low) feedback leads to a 
low (high) coronal temperature. This causes the shift of the spectral 
break to lower energies with decreasing /3. 

The moderate changes in the amplitude and direction of the 
bulk velocity thus induce crucial changes in the emitted spectra. 
Note that in most of our simulations, the bulk velocity is smaller 
than the thermal electron velocity /3th (e.g. /3th ~ 0.5 for kTe = 
100 keV). 



4.2 The ampliflcation factor 

The low feedback at /3 > iniplies a high Compton amplification 
factor of the X-ray source (eq. hi). The main effect can be estimated 
analytically assuming that the source is isotropic in the plasma rest 
frame (see eq. 7 in B99b), 
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Figure 6. Amplification factor A as a function of /9. The simulations are 
done for a cylinder with h/r = 2. The results are shown by symbols with 
the same meanings as in Fig. ^ The solid and dotted curves display the 
analytical formula for GBH models with tt = 0.5 and rr = 3, re- 
spectively. 



A{(3) = A,, 



' 1-/32(1 + ^,)V4 



(7) 



Here, 7 = - /32 and Us is the effective geometrical param- 

eter which we take from the static model (see eq. ^ and Fig. |3|b); a 
more exact definition of /is is given in B99b. A — Astatic at /3 = 
and monotonically increases with /3; at /3 < equation (Q) gives 
A<A 

static- 

In Fig. P we compare the resuhs of simulations with the an- 
alytical estimate (Q). In the optically thin models (tt = 0.5), we 
observe the reduction in A for large positive P compared to equa- 
tion (^). This is caused by the anisotropy of scattering in the plasma 
rest frame (optically thin models have very high temperature, see 
Fig. and scatter preferentially downwards). At /3 < 0, is 
small and the agreement is good. 

In the optically thick models (tt = 3), we see that the actual 
value of A at /3 < is smaller than that given by equation (Q). This 
is caused by the trapping of radiation and its advection downwards 
by the moving plasma. The advection enhances the anisotropy of 
the blob radiation compared to the optically thin case. In the lim- 
iting case of large optical depths, the velocity of photon diffusion 
PdiB ~ I/tt is smaller than the advection velocity, P, so that al- 
most no radiation can escape through the upper boundary. When 
viewed from the plasma rest frame, the effect can be understood 
by noting that both the bottom and top boundaries of the cylinder 
move upwards; as a result, the escape probability is larger for pho- 
tons emitted towards the bottom boundary. 

For the similar reason, equation (^ underestimates A for high- 
tt models at /3 > 0. In this case, the radiation is advected towards 
the upper boundary, and it is strongly beamed away from the disc, 
more than in the optically thin case. 



Figure 7. F versus A. The results of simulations for GBHs and AGN are 
shown by filled and open circles, respectively. The source geometry is fixed 
(a cylinder with h/r = 2), while tt and /3 vary. Only models with 50 < 
kTe < 300 keV are shown here. The lines correspond to approximation (^. 



4.3 The spectral index 

A simple functional shape T{A) suggested in B99a fits the results 
of the simulations. 



r = C(yl-l)" 



(8) 



The fitting parameters are C = 2.19, 5 = 2/15 for GBHs (fcTbb = 
150 eV) and C = 2.15, 5 = 1/14 for AGN (fcTbb = 5 eV). These 
parameters are not far from those found by B99a in the case of a 
spherical blob: C = 2.33 for both GBHs and AGN and 5=1/6 
and 1/10 for GBHs and AGN, respectively. Equation (^) is a good 
approximation for models with fcTe in the range between 50 keV 
andSOOkeV (Fig.^. 

Equation (g) combined with equation (^ yields the analytical 
dependence r(/3) which we find to be in good agreement with the 
results of simulations (see Fig. ^). Only for GBHs models with 
very high electron temperatures (kTe ^ 300 keV) the deviations 
AF exceed 0.1. The reasons of these deviations are the same as in 
the static case discussed in the end of Section 3.1. 



4.4 The amplitude of reflection 

Following the analysis of B99a,b, but accounting for the attenua- 
tion effect (see Section 3.2), we get an analytical formula for the 
reflection amplitude. 



R{p) 



+ 



(l-Mo)[l + /3(l + Mo)/2] 



(9) 



-tt{1-Mo) 



where, fi = cos i, and i is the disc inclination. Here, the reflected 
luminosity is represented as a sum of two parts: the first one is 
reflected outside the cylinder base and does not experience any at- 
tenuation and the second one is reflected from the base and it is 
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Figure 8. Characteristics of tlie dynamic disc-corona model. The model is determined by four parameters: h/r = 2, the optical depth of the cylinder, rx 
(measured in the vertical direction), the temperature of the soft radiation entering the blob, Tbb, and the velocity of bulk motion /3. The observed characteristics 
are inclination-averaged. Filled and open symbols correspond to fcT^b = 150 eV (GBHs) and fcT^b = 5 eV (AGN), respectively. Squares and triangles show 
the cases tt = 0.5 and rr = 3, respectively, (a) The average temperature in the blob versus /3. (b) The amplitude of reflection R versus /3. The curves 
give analytical /j-averaged R from equation (H) (dashed - T'v = 0.5, dot-dashed - rx = 3). (c) The photon index of the intrinsic X-ray spectrum F (in the 
2-10 keV band) versus /3. The analytical approximation (eq. ^ and ^ is shown by solid curves for GBHs and dotted curves for AGN. (d) Reflection R versus 
spectral slope F. The analytical model is shown by solid curves for GBHs and dotted curves for AGN (left curves for rx = 0.5 and right curves for rx = 3). 



partially attenuated, depending on rx- When /^o approaches unity 
(or rx — > 0), the attenuation is not important and equation (^ be- 
comes equation (3) in B99b: 



(l + /3/2)(l-/j^f 



(10) 



Equation (^ is in good agreement with the simulations (see 
Fig. ^). Substantial differences appear only in the case of low-rx 
GBH models. The differences are caused mainly by the anisotropy 



of scattering in the rest frame of the hot plasma. The similar effect 
was observed in the static GBH models (Fig. |^) and discussed in 
Section |3.4 

The left panel of Fig. ^ shows the dependence of reflection on 
the disc inclination. Note that R is less sensitive to P at large incli- 
nations. The overall angular dependence of R is quite well repre- 
sented by equation (^. 
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Figure 9. (a) Reflection, R, versus inclination (fj, = cosi). Squares are for r-r = 0.5, triangles for tt = 3. Here h/r = 2, fcT^b = 5 eV, and = 
—0.2, 0.1, 0.5. Reflection is larger that unity for negative f3 (bulk motion directed towards the reflector), and is much smaller than unity for the outflows. The 
dotted curves are the analytical estimates given by equation (bI). (b) R — T correlation for two inclinations /i = 0.25 (i = 76°, dashed curves), and fi = 0.85 
(i = 32°, solid curves). Here (3 varies from —0.3 (upper points) to 0.5 (lower points) with the step 0.1. The curves connecting the points are added for clarity. 



4.5 R — T correlation 

Fig. ^ shows R versus F. Again we see that the analytical formulae 
are in quite good agreement with the results of simulations. The 
R — T correlation is steeper in AGN than in GBHs. This is caused 
by the fact that r{A) is a flat function at small Tbb (see Section 
4.3; the index 5 in eq. ||is smaller for AGN than that for GBHs). 

The shape of the R — T correlation also depends on the disc 
inclination: the reflection amplitude is very sensitive to /3 at small 
inclinations and the correlation gets stretched in the vertical direc- 
tion (see Fig.|^). 

4.6 Pair dominated coronae 

A compact and energetic blob emitting hard X-rays may be mainly 
composed of pairs produced in 7 — 7 reactions. In this case, the 
blob optical depth and temperature are determined by the compact- 
ness parameter /diss ~ LdisBiyT / hnieC^ , where Ldiss is the power 
dissipated in the blob (see e.g. Stern et al. 1995a,b). The equilib- 
rium tt is controlled by the balance between pair production and 
annihilation and Te is controlled by the energy balance. 

We have calculated a set of pair-dominated dynamic mod- 
els with the h/r — 2 geometry. For simplification, we compute 
a global pair balance assuming a uniform pair distribution in the 
cylinder. The results are shown in Fig. 

At fixed l3, the model follows an equilibrium Te — tt curve. 
We emphasize that the Te — tt relation for plasma is the same 
as for normal e — p plasma, since it is determined by the energy 
balance only (see HM93; Stern et al. 1995b; PS96) and the annihi- 
lation radiation does not contribute much to the total energy budget. 

With increasing (3, the Te — tt curve moves to the right. An in- 
crease in 13 at fixed Zdiss leads to an increase in the optical depth 
while the temperature does not change much. This is the 'thermo- 
stat' effect of pairs discussed previously in the static models (see 
Figs 2 and 3 in Svensson 1996). The typical values kTe ~ 100 
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Figure 10. The relation between Te and -pr for /3 = -0.3,0.0,0.3,0.5 
(from left to right). The simulations are done for a cylinder with h/r = 2; 
the open and filled circles show the AGN and GBH models, respectively. 
The Te — Tx relation is shown by dotted curves for AGN and solid curves 
for GBHs. This relation does not depend on the nature of the scattering 
plasma. The dashed curves, desciibing the case of pure plasma, connect 
points of equal dissipation compactnesses idiss ■ 



keV and tt ~ 1 — 2 inferred from observations (Zdziarski et al. 
1997; Poutanen 1998) are obtained for large but possible values of 
the dissipation compactness, /diss > 100. 

We now briefly discuss the J? — F correlation expected at fixed 



© 0000 RAS, MNRAS 000, 000^000 



10 J. Malzac, A. M. Beloborodov and J. Poutanen 



/diss and varying (5. The correlation will be steeper compared to the 
tt = const case shown in Fig. |^d. At large /3, the equilibrium tt 
is large and here the R ~V track is close to that with tt = 3. At 
/? < 0, Tt drops and the R — T track approaches the tt = 0.5 
curve. 



5 COMPARISON TO OBSERVATIONS 

5.1 Broad-band spectra of individual sources 

Fig. [ll| shows the spectra of one GBH source, Cyg X-1, and one 
bright Seyfert 1 galaxy, IC 4329A, together with model spectra. 
We did not use a real fitting procedure (since each simulation takes 
a few hours on a modern workstation), instead, we used unfolded 
spectra and searched for similar spectra in our sample of models. 
We find that the spectrum of Cyg X-1 is well reproduced by the 
model with a bulk velocity /3 = 0.3, an optical depth tt = 3, 
and h/r = 1.25. The inclination that gives the best agreement 
with the data, i — 50°, is compatible with the current estimates 
25° < i < 67° (see Gierlinski et al. 1999 and references therein). 

For IC 4329A, we find /3 = 0.1, tt = 3, and h/r = 2. 
The inclination, i = 40°, is consistent with the Seyfert unification 
scheme stating that Seyfert Is should be seen rather face-on (An- 
tonucci 1993). Note that h/r is smaller in the case of Cyg X-1. 



5.2 R — T correlation 

Fig. |l^ gives the observed R — F correlation for Seyferts and GBHs 
(Zdziarski et al. 1999; Gilfanov et al. 1999). The correlation exists 
in individual objects observed at different periods as well as in a 
sample of sources. Such a correlation supports the hypothesis that 
it is the reprocessed/reflected radiation that cools the X-ray source 
via Comptonization. 

The R—T correlation expected in the dynamic corona model 
is similar in shape to the observed one, and we now try to fit the data 
with the model. As we showed in Section 4, the analytical model is 
a sufficiently good approximation to the results of exact simulations 
and it is reasonable to use this simple model in a fitting procedure 
instead of the exact transfer model (especially taking into account 
the uncertainties in the validity of assumptions such as homoge- 
neous density and uniform heating of the blob). 

The analytical model yields explicit expressions for R{f3) 
(eq. ^ and r(/3) (eqs. ^ and |^ for given parameters /is, Mo, a, 
and Tt. For simplicity, we assume /j.^ = Mo and fix the albedo at a 
reasonable value a = 0.15. The fit to i? — F correlation is shown 
in Fig. |l2|for AGN and GBHs. 

It is interesting to note that the correlations for AGN and 
GBHs cannot be fitted with the same value of /is. We find fis be- 
ing lower for GBHs (/is ~ 0.45 or equivalently h/r — 1) than 
for AGN (/is = 0.6 or h/r — 1.5). This might indicate different 
geometries of emission regions in GBHs and AGN. 



5.3 Tlie beaming effects 

The observed flux Fh of the hard X-rays beamed away from the 
disc can dominate over the blackbody component _Fbb when the 
disc is observed at small inclinations. The hard X-rays can dom- 
inate even if only a modest fraction / of the accretion energy is 
dissipated in the corona. In the case of a patchy corona, the model 
predicts 



i?(i)2cosi 1 



1-f 1+d 



f 



d 



(11) 



Here d — L-/L+ is the ratio of the Comptonized luminosities 
emitted towards and away from the disc. 



d: 



1 + /3/2 (1-/3)= 



-/3/2 (l + /3)= 



(12) 



For example, in the case of Cyg X-1, taking i = 50° and f3 = 0.3, 
one gets F^b/Fh ~ 0.4 for / = 3/4 or Fbb/Fh « 0.85 for 
/ = 1/2, in agreement with the data (e.g. Di Salvo et al. 2000). 

The beaming of the coronal X-rays may be tested in a sample 
of hard-state objects with known inclinations: R should increase 
with i. Also, the shape of i? — F correlation varies with i (Fig. |^, 
which may be checked in the data on GBHs. Unfortunately, in the 
case of AGN, the high inclination objects are obscured by the large- 
scale molecular torus which dominates the reflection component. 

The ejection model with [3 < predicts very soft spectra and 
large ratios -Fbh/^ii- This situation might take place in some soft 
objects, e.g. Narrow Line Seyfert 1 galaxies (see Brandt 1999 for a 
review). 



6 DISCUSSION 
6.1 Variability 

The coronal emission should be variable, with a temporal pat- 
tern governed by the complicated processes of the magnetic field 
buoyancy and its dissipation in the corona. The ejection velocity 
can vary from flare to flare and the fit parameters of the observed 
spectrum should be taken with caution since they give only effec- 
tive/average values. 

The main X-ray emitting region probably extends through a 
decade in radius, and l3 is likely to be different at different radii. 
Recent temporal studies (e.g. Gilfanov et al. 2000; Revnivtsev et 
al. 2000) suggest that harder spectra with lower reflection come 
from the inner radii. In the context of the dynamic corona model it 
would correspond to an increase of [3 in the innermost region. The 
analysis of Cyg X-1 and GX 339-4 at different periods of the hard 
state showed that the typical variability time-scales grow when the 
spectrum hardens. This behaviour may indicate that the accretion 
is slower in the hardest periods (Beloborodov 2001). 



6.2 Pair envelopes in quasars 

In this paper, we neglected scattering on the particles escaped or 
created in 7 — 7 collisions outside the heated region. In the bright- 
est sources, with luminosities comparable to the Eddington limit, 
the hard X-rays from the flares will interact with each other and 
create an optically thick outflow covering the whole inner re- 
gion of the accretion disc (Beloborodov 1999c). Such an outflow is 
relatively cold (at the Compton temperature, kTc ~ 1 — 10 keV) 
and it has an equilibrium bulk velocity such that the radiative pres- 
sure from the disc is balanced by the Compton drag. The velocity 
is mildly relativistic (l3 increases from 0.3 at the base to 0.7 at the 
photosphere of the outflow) and the scattering of the disc/corona 
radiation in the outflow strongly affects the angular distribution of 
the observed luminosity. The reflection amplitude will then be sup- 
pressed whatever dynamics takes place in the compact flares hidden 
in the outflow. 



© 0000 RAS, MNRAS 000, 000-000 



X-ray spectra of dynamic coronae 1 1 




E (keV) E (keV) 

Figure 11. (a) Spectrum of the Seyfert 1 galaxy IC 4329A observed by ROSAT, Ginga, and CGRO/OSSE (crosses, data from Madejski et al. 1995). The low 
energy cut-off is due to the Galactic absorption. The solid curve shows the model spectrum for rr = 3, h/r = 2, /3 = 0.1 at inclination i = 40°. (b) 
Spectrum of Cyg X-1 as observed by Ginga and CGRO/OSSE in 1991 September (crosses, set 2 from Gierlinski et al. 1997). The solid curve shows the model 
spectrum for tt = 3, h/r = 1.25, /3 = 0.3 at inclination i = 50°. In both panels, the dotted curves give the reflected components and the dashed curves 
show the intrinisic Comptonized spectra. 
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Figure 12. (a) Thepl^erved R — V correlation for Seyfert galaxies (circles). Data from Zdziarski et al. (1999). The curves show the correlation predicted by 
the model (eqs. pi Wfl). Here /Xs = /lo = 0.6, r-p = 3, a = 0.15. The solid, dashed, and dotted curves correspond \o fi = 0.85, 0.55, and 0.35 (i.e. i = 32°, 
57°, and 70°), respectively, (b) The observed R-T correlation for GBHs, Cyg X-1 (circles), GX 339-4 (triangles), and GS 1354-644 (squares). Data from 
Gilfanov et al. (2000). The curves correspond to the same three inclinations as on the left panel. Here fia = Ho = 0.45, tx = 3, a = 0.15. 



6.3 Optical and UV polarization of Seyferts and quasars 

The scattering of soft radiation in a dynamic corona is accompanied 
by a strong effect that can be tested in optical/UV observations: the 
polarization of scattered radiation is sensitive to /3. The polarization 
is perpendicular to the disc if 0.12 < (3 < 0.78 and parallel oth- 
erwise (Beloborodov 1998). The observed optical polarization in 
non-blazar AGN is parallel to the radio jet that is presumably per- 



pendicular to the accretion disc (Stockman, Moore & Angel 1984); 
it is consistent with the scattering in a mildly relativistic outflow. 

In three of ten quasars with measured UV polarization, a steep 
rise in polarization was observed blueward of the Lyman limit (see 
Koratkar & Blaes 1999 for a review). The outflowing corona model 
may provide a natural explanation for this mysterious rise (Be- 
loborodov & Poutanen 1999). 
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7 CONCLUSIONS 

We performed detailed Monte-Carlo simulations of X-ray produc- 
tion by hot coronae atop accretion discs, and tested the model 
against the data. The main results are as follows. 

(i) A static corona atop a neutral reflector is not able to pro- 
duce the observed hard spectra with low reflection. Furthermore, 
changes in the coronal geometry produce an anticorrelation be- 
tween the hard X-ray spectral slope F and the amplitude of reflec- 
tion R. This anticorrelation is opposite to what is observed. 

(ii) The disc-corona model becomes consistent with the data if 
the hot plasma is dynamic, i.e. moves with a mildly relativistic bulk 
velocity, 13 — v/c, with respect to the accretion disc. In particular, 
the spectrum of Cyg X-1 is reproduced by the model with /? = 0.3, 
confirming the estimate of B99b. 

(iii) The observed R—V correlation is well explained by varying 
/3. It suggests that P may be the main parameter controlling the X- 
ray spectrum. 

(iv) The results of the simulations are in good agreement with 
the analytical description of B99a,b. We improved the analytical 
model by accounting for the attenuation of the reflection compo- 
nent by the hot plasma atop the disc (eq. ^). 

(v) The optical and UV polarization properties of AGN are sen- 
sitive to the dynamics of the corona. The polarization data provide 
an independent indication for the mildly relativistic bulk motion. 
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